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ABSTRACT 


The advent of tilt rotor technology asks for rotors that 
have different twist and RPM requirements in hover and in 
forward flight to optimize for operational conditions. In 
order to get an assessment of the capabilities to fulfill 
these requirements this report presents a mapping of twist 
angle variation as a function of RPM and laminate orientation. 
The basic laminate for the six models as well the D-shape spar 
that represents the structurally active part of the blade is 
assumed to be constant (0°/90°/0°/90°/90°/0°). This six laver 
cross ply laminate is chosen as it provides the necessary 
extension-twist coupling without a hygrothermally induced 
twist that is highly undesirable. The couplings and trends in 
the models are visualized in carpet plots, one for each model, 
in an attempt to establish a method to answer the basic 
question of the magnitude of twist angle available due to a 


particular geometry, material and load system. 
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[ CNTRODUCTLION 


A . OVERV IEW 


In the design process of composite material structures it 
is often noticed that the best characteristics of the material 
ere not used. The normal procedure is to tailor the material 
to approach the properties of isotropic materials, avoiding 
the ‘undesirable’ couplings between extension and twist for 
example. 

That traditional approach worked for helicopter design. 
However, in the tilt rotor, the idea of changing the twist 
fee Of the blade to optimize for either hover or forward 
flight came into the scenario. 

Mors Stwa. will wse a blade spar to study the sensitivity 
Ore cthne extension-twnist coupling due to changes in RPM, laver 
angles definition, changes in laminate orientation and also 
the effects of each of these configurations in the first four 
natural vibration modes of the blade. 

Mere allaiIxciS, a program named CASA/GIFTS (Compuler 
Pemeec otrictural Analvsis/Graphical Interactive Finite Element 
Meee oxetem} is adopted. The program makes it easier to 
Memeiaeerae the effects of the coupling in the structure and is 
Memeo fe als C€lficient design tool. The resulis obtasmed 


Mreousl: the Moe Lemetit procedurt ee. each laminite 


configuration, in a total of six models will be presented in 
3-D graphs that show the trends for each model for the laver 
angle (top layer as a reference), the RPM and the twist angle. 

These graphs can then be used as a preliminary orientation 
by the designer to achieve the optimal or near optimal 
combination of laminate configuration layer angle within the 
laminate that fit the specific requirements for the particular 


blade in focus. 


de i LITERATURE REV ITILEW 


This section reviews briefly some pertinent publications 
listed as References in the thesis report. Other related 
literature found during the scope of this research, is listed 
under Bibliography. The topics reviewed are classified under 
the following categories: 

Finite Element Theor. 
Composite Material Theory. 
None a rons. 

Helicopter Theory. 

Design Applications (Modeling). 
Auxiliary Software. 

In the Finite Element Theory, relevant to present 
research, some basic references are: Batoz[1980] thane 
describes the finite element procedures in engineering 


analvsis using Discrete Kirchoff Theory (DKT), Cook/ 193 asian 


basic concepts and applications of finite element analvsis, 
Craig[1981] describing computer methods in structural dynamics 
and [CASA/GIFTS, 1987] that is the reference manual for the 
finite element code used for the present (static and dynamic) 
work. Its graphical capabilities are very convenient to the 
analyst ina manner that simplifies the visualization of the 
effects of loads and couplings in the structure, making 
difficult concepts more accessible to understanding. 

Tsai and Pagano[1968] establish a notation in which the 
composite lamina properties are invariant with respect to the 
axis of rotation. This approach provides a verv useful way to 
compare various material systems; the laminate theory is well 
documented, see for example Vinson[1987] and Jones[1975] both 
offer a good source on the mechanics and behavior of composite 
materials. 

Yntemal1955] developed a very useful tool to estimate 
bending frequencies of rotating beams in which the stiffness 
effect of RPM is shown for several blades (beam) 
configurations. Schilhansl[1958] uses another approach to 
estimate bending frequencies of rotating cantilever beams, 
introducing the effect of the angle made by the minor axis of 
inertia with the direction of the tangential circular 
velocity. Wood[1965]|] shows a parametric investigation of 
aerodynamic and aeroelastic characteristics of articulated and 
hingeless rotor systems. Pritchard][1988] presents an optimal 


placement of tuning masses to reduce vibration levels in 


helicopter blades. Kottapalli{1983] approaches the vibration 
reduction problem by modifving the blade torsional response. 
This 1S an innovative approach in that he modifies the 
aerodynamic response instead of adding mass, which is the 
usual conventional way to solve the vibration reduction 
problem. 

The helicopter theory in general is found comprehensively 
in Johnson{1980}] and is complemented in this research by 
Prouty[1986] in the understanding design trends. 

Hoskin[1986] is a general informative book in the 
aeronautical utilization of composite materials, and in 
Lake[1988] there is a preliminary investigation of finite 
element modeling Of composite material rotor blades. 
Nixon[ 1987 } addresses the extension-twist coupling oy 
composite circular tubes applied to tilt rotor design. In 
McVeigh [1983] the aerodynamic design of XV-15 advanced 
composite tilt rotor blade requirements are reviewed. Finally 
Hodges/]1987] presents a comparison of composite material rotor 
blade models, using two different methods of analysis: this 
model was chosen as a prime reference to build the basic model 
used throughout this work. 

Several software tools made it possible to manipulate the 
immense volume of data transforming, it into a suitable and 
compact presentation form, enabling the conclusions to be 
obtained. kKeliv{1988] was the basic reference to use Word 


Perfect 5.0 in the word processing and graphical! editing sees 


Holt{1988] is a reference that helps to manage the data into 
Lotus 1-2-3 files and to get the input graph files to Word 
Perfect and to SURFER[1987]. The SURFER program was used to 
obtain the 3-D carpet plots of extension-twist coupling 


analvsis. 


— . THESIS OUTLINE 


This thesis report is divided into three areas, often 
found during a design process. 

It begins by modeling a rotor blade with laminate 
composite material. A numerical solution i1is-~ obtained, 
exploiting the extension-twist coupling inherent when 
asymmetries are present in the laminate construction. 
Discrepancies encountered in the results compared with other 
rotor models are reported. These differences raised questions 
fostering the research towards the possibilities of different 
rotor constructions, starting with the same basic laminate. 

The different rotor configurations obtained are then 
subject to both static and dynamic analysis. The static 
analysis reveals the different responses of these models, 
while the dynamic behavior remains almost invariant, vielding 
only small variations in the first four fundamental vibration 
modes. 

Analvsis of the results is done, and it is believed that 


a new insight into finite element analysis of composite 





laminates is provided. This new insight is bett 
there are asymmetries present within 


construction. 


1 ae THEORETICAL FORMULATION 


Ae BAS TC CONS TDERATIONS 


To understand how composite material works and its 
particular behavior that is useful in "designing" the 
material, a macroscopic approach is resorted to. 

Composite means that two or more materials are combined 
and usually preserve their best qualities, and, often exhibit 
properties that neither one possesses by itself. 

Some of the properties that can be improved by making a 


composite material are: 


Strength 
stiffness 
corrosion resistance 
. wear resistance 
. life in fatigue 
». thermal insulation 
thermal conductivity 


acoustical insulation 


It may be noted that not all these properties are improved 


at the same time, and usually it may not be required. 


~] 


An Orthotropic! material has properties which are 
different in three mutually perpendicular directions at a 
point, in other words has three different planes of symmetry. 

A lamina is the basic form of any laminated fiber- 
reinforced composite, it consists of a plane mounting of 
fibers imbedded in a matrix. 

In the macro-mechanical behavior of the lamina, only 
averaged apparent mechanical properties, or better, the 
stress-strain relations for an orthotropic material under 
plane stress conditions, transformed to directions not aligned 
with the principal directions of the lamina, are considered. 

This describes the behavior of the orthotropic material 
whose laminae have arbitrary directions, other than the 


natural geometric directions of the structure to be analyzed. 


1S STRESS -STRAIN RELATIOONS 


The stress-strain relations or generalized Hooke’s law for 
anisotropic material will be reviewed briefly, and the basic 
restriction of linear elastic behavior assumption applies. 

The following derivations can be found in more detail in 
Vinson! 196 7c. 

The elasticity tensor has nine independent elasticity 


constants, because of symmetry (C,,;=C;;) and the three mutually 


1_- Orthogonally anisotropic. 


8 


orthogonal planes of symmetry (C,g=C2,=C3g=C45=0 ) ; 


what is, 
C13 Cy2 Cy3 9 0 0 
C21 C22 C23 9 O DO 
ei) oA; = C3) C32 C33 0 0) 0 
0 Cag O O 
0 0 0 O Ce. O 
| 0 0 0 0 O Cee | 
and 
(2) Cree Wc ee lee, 6}) j=1,2,3,4,5,6)* 


In these relations the hygrothermal and thermal effects 
are not included for the sake of simplicity of this report. 
Through a force equilibrium study, analogous to the Mohr’s 


circle analvsis [Jones, 1975] the stresses and strains in the 


‘The material directions 1 to 6 are equivalent to 
cartesian directions: x*,Vv,Z,VZ,XZ,xv respectively in material 
axis. 


material directions can be related to arbitrary’ directions 


uSing the following relations: 


Oy O71 

Oo oO 

y 2 

(3) oO, =( Te alos 
Oyz O4 

Oxz Os 


| Oxy [6 


and 
ex Ey 
ey E2 
(4) E, = ile Ee, 
Eyz E4 
Exe E5 
Exy | re 


The transformation matrix is constructed by defining 
m=cos 86, n=sin 8; where 8 is the ply orientation angle, 
defined with respect to material and structural axes, assuming 


positive sign in counterclockwise direction (Figure 1). 





3 More convenient, xyz axis related to the structure 


geometry for example. 


10 





Figure 1 Positive rotation of principal material axes from 


arbitrary xy axes [Jones,1975]. 


The inverse of transformation matrix [T] is expressed in 


terms of direction cosines m and _n as: 


ms ns 0 0 0 -2mn 
n m 0 0 0 +2mn 
-1_ 6) 0 ) 0 0 0 
(5) [T} “= 0 0 0 m <-n 0 
0 0 0 n m 20 9 
-mn 0 0 0 0 (m - n ) 


tal 


For a lamina of composite material the modern notation 
uses Q;; for the quantities in the material stiffness matrix 


and are defined as: 


QFE (1-Vo3Vg9)/8) = QagF Egg (1-V31V,3)/6 


)/5, Q = G Q G 


44> 793? 55 “ie 
(6) QFE 6% 91% 31% 93 )79 = Eon l%y 04% 30% 3)/8 
O37 Ey 1 6% 314% 01% 30)75 = Eno l¥y3t¥ 10% 93)/8 
Q 537 Eo (Vg0t¥y0%31)/5 = EgglVogtVo V4 3)/5 

6 = 1 — Vi9%917 Y23%307 Y31%137 2%21%32%13 


If the lamina has same properties in both 2 and 3 
directions, then vy2=Vy3)  V3)=V32=0, G,2=G,3, E22=E3, and a 


simplification in equations (6) lead to equations (7). 


Q54 FE y/'l - VyoVo,): QooEgo/(1 - Vy 9Vo1! 
(se) Q,07 Vo E,,/1 - Vo. ont Vi9Eo9 /(1 - V49%}1? 
667 S12 


iia 


To compare various material systems in the design of a 
composite structure, Tsai and Pagano[1969] arrived at some 


invariants with respect to the axis of rotation: 


cee 
Uy= g (3Q,,+3Q55+2Q, ot4Q¢,) 
U.= = (e .- ome 
o= 3 (811- Poo 
ae : a 
(8) Us= § (Q),+ Qoo- 2Q,5- 4Qee) 
Ue Quet Gr Gone) 
a= g §9117 200 12 66 
U-= = (Q..+ @,,- 20,.+ 40...) 
5= g (844+ Qo 12 66 


These invariants are representative of the lamina 
properties. The [Q] matrix may be transformed to x-y axes from 


1-2 axes by the usual coordinate transformation. 


1 


(9) fale = it) mrontn | 


In terms of the invariants Us, the [Q] matrix elements 


are written as: 


aa U,+ U,cos( 20) + U,cos(46@) 
Qo0> U,- Up,cos( 26) * U,cos(4@) 
Q,07 U,- U,cos(46) 

ier ) - | 
Qe6= U; U,cos(46) 
e _ 1 ; 
Qi¢ 5 U,sin( 26) + U,sin(40) 
Q.-2 ‘ U,sin(260) - U,sin(40) 

26 Z 2 ys S 


als 


After the definitions in (10), the stress-strain relations 


for the kt® lamina of aN ply laminate (Figure 2) are given by: 


a ex 
"y oy 
= = 
(11) a2 | = [@l, | e 
ous ee 
eee a4 
XY J, L *Y J, 


Where @,,; are defined earlier. 
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Figure 2 Laminate construction, [90°/+45°/0°/-45°/90° J 
characterized in this case to show the nomenclature for 


the lay-up. 


ie 


om LAMINATE ANALYSIS 


When two or more laminae are bonded together to obtain an 
integral structural component, the laminae principal material 
directions are established to produce a component that can 
resist loads in several directions (Figure 2). 

The following procedures enable to obtain the stiffness 
matrix of such a composite material from the lamina 
properties, using arbitrary combinations of materials and 
layer angles. 

For the kt® lamina of a N layer laminate, equation (11) 
represents the stress-strain relations. 

In order to obtain the strain-displacement relations, the 


displacements are defined as follows: 


i a Ga oye) 


Ulex ay 2 


(12) Vil oe eee) ey 2 ( x) 


Bul Zz) w(x,y) 


where u, v andw are displacements in x, y and z directions 


respectively 


LS, 


and? 


5w bw 
(13) qa = =——— and 3 = -— 


5x Sy 


are the negative of the slope, related to the rotations 


— 


The mid-surface strains, then are given by 








Ee Su, “ E.. ON ’ 
=o OBS “9 7 by 
(14) 
Ee  ._ 1 Su bv 
SO Oo.) UF fe) 
ony 5x 


while the curvatures are given by 


Yu In the following expressions 6 means pamt jal 
derivative. 
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Defining h as laminate thickness, h, (Figure 3) as the 
vectorial distance from the panel mid-plane, N, M, Q@as stress 
resultants (Figure 4), stress couples and shear resultants 
respectively. Relations can be established relating these 
resultants to the strains and curvatures. This matrix is 
henceforth referred as the ABD matrix. 

This nomenclature is adopted for clarity and simplicity 
as may be noted from the expressions that follow. 

The stress resultants are expressed bv integrating 


stresses across the thickness of the plate and is given by: 


N e) 
% x 
Me +h/2 oy 
(16) Ne = | oy dz 
Q . -h/2 a 
x XZ 
Q., | os 


and the moment resultants take the following form 


x +h/2 Ox 
el 7 } M = | Oo Z az 
-h/2 v 
Oo 
res aay 


ae 





Figure 3 Nomenclature for the Stacking Sequence 


[Vinson,1987]. 





Figure 4 Positive directions for Stress Resultants and 


Stress Couples for a Plate [Vinson,1987]. 
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Relating resultants with the inplane strains” and 
curvatures through expressions (11) through (14), matrix form 


1S given by: 


[N ] fone) + Bae | K] 


(18) 
[IM] = [B].[eE,} + (D).[K] 
In expression (18), the matrices A, B and D are called 
extensional stiffness matrix, bending-stretching coupling 


matrix and flexural stiffness matrix respectively, where: 


N raat s e 
- ¥ _ = 
1 N = 2 2 ot 
il N : 27... grms a 
Des oe * ea a 4 Ne ny ees [i,J=1,2,6] 


Assuming that transverse shear stresses are distributed 


parabolically (even though there exists discontinuities at 


the laminae interfaces), the shear resultants may be written 
as: 

@.. = 2( Az zEL, + Sar 
fez 0 } 

Q. = 2( Agse, - news 


where 


a 9 e oe Se 1 
2a) A oe ( See hi [ he + aa )* "| 


1,jJ = 4,5 only 
The stresses and moments resultants can be written in a 
compact form that is very convenient to see the overall 
relations within a composite material laminate, specifically 
the ABD matrix. The ABD matrix is obtained as a combination 
of the expressions in (18). Expressions (18) and (22) are 
equivalent and in expression (22) the ABD matrix is defined 


explicit lw 





Ne Ebest Et AQ 2 _ 

Ns per we 2 Pim Bao - “a 
(22) Nxy | = | 416 426 See : Pie Poe Peg | * | Sxy, 

MS. Bia Pio Pile ey ene K. 

MY Bil 22 "260 eenie eee KY 

Mey Bis Bog Beg : Pig Pas P Ky 
1m EFFECTS OF PLY ANGLE e 


To get an insight of the overall influence of the ply 
angle 6 on the laminate stiffness, it is useful to observe 
the effects of 6 for a given lamina. As each element of the 


B. 


ABD matrix, A ij, and D,;, is a function of the material 


ij? 


stiffness coefficient Q: 


i;: the layer thickness t, and mid-plane 


ho 
S 


vectorial distance h,, the Q;,'s represent the stiffness 
variations due to ply angle changes ina given lamina. 

The property of the composite material that reveals its 
advantages over metallic or plastic is that it can be tailored 
to fit a given structural shape and a specific set of loads 
and operating conditions. 

The terms of the ABD matrix are functions of the material 
properties and the ply angle 0. If the properties are kept 
constant and within typical values of most fiber resin types 
of composite, it is possible to plot the values of Q,,; versus 
6 and get a set of seven verv useful graphs. These will permit 
an initial appreciation of such effects, and how they can be 
combined to meet a specific structural requirement. 

Figures 5 to 11 show how the elements of the "material 
stiffness" matrix ABD can be tailored by changing the ply 
angle, taking full advantage of this characteristic behavior 


of laminated composite material. 
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In figure 11 all Q’s are superimposed, on the same plot, 
to give a better perspective of the trends in doing this kind 
of laver design to fit specific structural constraints. When 
Aig and Agg are non-zero there is stretching-shearing coupling, 
while non-zero B,g means twisting-stretching coupling, and non- 
zero Bog indicates bending-shearing coupling; bending-twisting 
coupling is originated in non-zero values of Djg and Dgg. 

In normal design procedures, these terms are made zeroes 
by selecting appropriate stacking sequence; cross ply (0°/90° 
plies combination) or angle ply laminates (-@0/+8 plies 
Somoynation). Other important factor to be considered is the 


Symmetry with respect to the mid-surface plane. 


The stress strain relations or Constitutive relations 
matrix), combined with the proper equations of equilibrium and 
the strain-displacement relations form the basis for analyzing 
beams, plates or shells. 

These basic relations lead to the motivation and reasoning 


behind the model tailoring to fit a design requirement. 
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nT ee ROTOR MODELING 


A e TYPICAIL CURRENT ROTOR 


MODELING 


Current rotor design is a compromise between the different 
requirements in hover and forward flight, and the analytical 
models reflect that tendency or limitation. For instance, the 
optimum twist for hover is not the same in forward flight for 
tilt rotors as well as for helicopters. Indeed, very often, 
the designers avoid extension-twist coupling with laminated 
compoSite materials because of certain undesirable effects. 
The result is that one of the two flight modes will have less 
than optimal propulsive efficiency, depending on design phase 
decisions driven by requirements and/or aircraft mission 
profiles. 

With the tilt rotor concept of mixing fixed wing and 
moLorcraft technology, these design trends were brought to 


surface once more. 


B. PROPOSED SOLUTION 


Nixon[{1987] proposed a solution for this sort of dilemma. 


The model used an extension-twist coupled rotor blade with 


aa 


100% RPM? in hover and 80% RPM in forward flight. These are 
reasonable ratios between hover and forward flight for a tilt 
rotor type aircraft. The difference in RPM provides a change 
in centrifugal force, inducing a change in the twist angle 
distribution, which can be used to optimize the rotor for both 
regimes. 

The use of composite materials makes it possible to design 
a rotor blade that takes advantage of the changing loads, 
rendering the required twist distribution optimum in each 


flight mode. At this point one question arises: 


What is the magnitude of twist available from a particular 


geometry and material system? 


The answer to that question lies in a method of 
prediction, with a model that behaves as close to the actual 
blade as possible vet inexpensive enough to permit a wide 
Spectrum parametric study, providing a source to evaluate 
several geometry and material systems. Certainly such modeling 
should be verified by experimental data. 

This research is aimed to address the analytical part of 
the above question and provide design tools to rotor blade 


systems. 


> Reference value for 100% RPM is 217 rad/s, in that 
sSscule of rotors Dlade- 
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The torsional moments that affect a rotating blade can 
modify substantially the twist distribution. 

In this numeric model, the pitching moments generated by 
aerodynamic loading are not taken into account, the emphasis 
being on the inertia and extension-twist coupling as sources 
of torsional moments. 

An assumed lift distribution acts on the model, but 
compared with the inertia forces, its effects are negligible 
on the rotating blade. The inertia forces acting on the 
rotating blade produce extension which causes twist by 
structural coupling. However, the same inertia forces produce 
a centrifugal flattening [Prouty, 1986] effect called tennis 
racket effect (because of the tendency of a tennis racket to 
align its plane with the plane of rotation as it is swung in 
em arc). 

The forces acting in a helicopter blade can be seen in 
Figure 12, and it may be noted that in this case the effect 


is reduced with balancing weights (Chinese weights). 


C. Seah CAPABITLETLES 


Gifts capabilities include generation of any type of 
structural model, with the following elements: 
one-dimension rods or beams 


two-dimension plates 


ag 


. three-dimension shells, solid and complex stiffened 

shells. 

The elements can be selected from a library of options. 

The materials can be created by the user or selected from 
a library of standard definitions. The program allows the user 
to define isotropic materials and laminated composite 
materials, that can have up to a maximum of hundred lavers. 

Loads can be concentrated or distributed, static as well 
as function of time. 

Static and dynamic analysis are performed on the models 
generated. The dynamic analysis provides free vibration 
response and forced response (transient computations), and 
model superposition. 

In using GIFTS as the finite element program to model the 
blade, the effects of the centrifugal force are taken into 
account with the forces acting on the lumped masses as shown 
in Figure 12, radially from the center of rotation. In using 
other finite element programs this feature must be verified 
for correct usage and application. The tennis racket effect 
can affect the twist angle to a large extent. If the 
centrifugal forces act in parallel lines instead of passing 
through the center of rotation, there will be no centrifugal 
flattening effect and results in larger twist angles. 

The basic rotor model selected to apply the proposed 
method is the one analyzed by Hodges/[1987]1, this model was 


analvzed using two different methods, namely by finite element 
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(MSC/NASTRAN ) and a 


coupled-beam model 
oe enter of Actation {[Rehfield, 1985]. The 
rot or geometry 
Blade Feathering Axis 
Chinese 
eel Presented by 
Centrifugal Hodges[1987 ]} was 


Force Acting on 
— Chinese Weight 


modeled in GIFTS using 


triangular TB3 element 


Centritugal (" flat triangular 
entniuga 
Force Acting 
Perse plate ee SE 
men 
element, including 


membrane, out of plane 


and Frctici16cus in- 
plane bending 





End View stiffness.’ [Batoz, 
1980]). The results 
indicated a wide 


Figure 12 Twisting Moments Due to disparity of results. 
Centrifugal Forces [Prouty,1986]. However, the present 
model appears to be 
within the assumption of linear and small displacements 
theory. 
The models in Figure 13 have a notation (M120 / M1-20) 
that is described in the next section. A twist angle 
distribution discrepancy could be explained by the tennis 


racket effect probably not present in the model described in 


CBl 


Hodges[1987], but the tip vertical displacement due to lift 
and weight with no rotation could not be explained by this 
argument. It may be observed that the displacements quoted in 
the referenced work are about 25% of the length, which appears 


to violate the linear theory that was used. 
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Figure 13 Comparison between GIFTS model and NASA TM 89024 


results. 


This discrepancy mav be due to a possibility that the 
models compared are not exactly the same, although every 


effort was made to take into account all the pertinent 


By 











geometric and material data reported. However this discrepancy 
was the driving factor that raised some questions, such as: 

What are the possibilities to model different structural 

configurations starting with the same number of layers, 

same lay-up design, same thickness and material 

Properties? (perhaps by inadvertent oversight by the 

designer, analyst using canned programs, or manufacturing 

engineer. ) 

The answer to this was based on the following assumptions: 

». The structure is a D shape spar. 

It is built using three shells, top, bottom and rear 

Shell surfaces 

The top and bottom meet at the leading edge and the rear 

Shell closes the box. 

PuGemeevaotished, the internal lay-up in the stacking 

Sequence remains constant. 

Prtne those assumptions a laminate to constitute offe shell 
meeece Nias be defined, using only the Laminate Principal 
MPementation (LPO) angle. For example, in the following 
mmmenate | +20°/-70°/+20°/-70°/-70°/+20°|, the top laver angle 
Memeoenoted +20° and is a six laver asvmmetric laminate. 

Rotating the whole laminate =O clockwise gives 
Meee 7 + 7 0°/-20°/+70°/4+70°/-20°], whose Laminate Praneipal 
Orientation angle is denoted as -20°. 

This laminate has intentionally designed asymmetr: to get 


Meee CeTISION-twisi coupling. Associated *ith this laminate 


is a normal vector whose orientation is seen in Figure 14. The 
D-spar then is built using three such laminates, with each 
laminate mounted with the normal vector facing inward or 
outward. This method of construction results in six different 
configurations of the rotor. These configurations and other 
variations within each one will be developed in the next 


section. 


D. VARIOUS ROTOR DES ITGNS OR 


CONE IT GURATIOONS 


The six different configurations are designated Ml to M6 
and the LPO®, varying from -90° to +90° in steps of 20°, 
provides the suffix to complete the model designation M1-90’, 
M340, M6-60 etc. In a total of 11 LPO in each model. 

Different configurations may be visualized with the help 
of the scheme shown in Figure 14. 

Model I and IV can be built with only one laminate wrapped 
around, one inverse of the other. 

In the program GIFTS, the outward vector for each surface 
is defined during the "GRID" definition, e.g. ina four sided 


grid, each shell is defined within a grid. The local xX axis 


6 Laminate Principal Orientation Angle, referred to the 
top laver of the laminate used to construct that particular 
model. 


7 M1-90 meaning model 1 with LPO = -90°; M340 meaning 
model 3 with LPO = 40°: 
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is oriented parallel to the first line of the grid, the local 
Z axis is the cross product of the first and second "line" in 
the grid definition, while the local Y axis is the cross 
product of the Z and X axes, as is usual in cartesian systems. 

The outward vector is oriented in a direction defined by 
the local Z axis. Layer angles are positive in the counter- 
clockwise direction in the XY plane for shell elements, as is 


the convention for laminated composite materials. 
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LAMINATE STACKING SEQUENCE 





Figure 14 Normal to surface vector associated with the 
layer sequence of the basic laminate and its position in 


relation to end view of the blade. 


It may be observed that, when an asymmetry is present in 


the basic laminate (material definition), the structural 
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response is highly dependent on how the grid and consequently 
the lay-up is defined. In the case of a symmetric laminate or 
isotropic material, this fact is never brought to surface 
because of its irrelevancy. However, dealing with asymmetric 
laminated composite materials, the analyst must always have 
this in mind. It may be noted that, in most cases the laminate 
1s chosen to be symmetric. 

Table 1 gives ina compact form the description of various 
laminate constructions described. Included in that table are 
the line names actually employed in the input file M$**.SRC 
[Annex 1] ($ = 1,2,3,4,5,6 models; ** = top layer angie 
defined previously, in the range -90 to +90), to develop the 
models generated in GIFTS. By modifying two areas in this 
input file, six models are generated that are considered under 
the present study, this also provides variation of the laver 
angle within each model. The actual modifications employed in 
the input file are given in LAYUP.AUX [Annex 2]. For each 
model (one of the six combinations obtainable by the change 
in the laminate associated vector, as given in Figure 14), 
the whole laminate was rotated in increments of 10°, 
referenced to the top layer angle (LPO). For each of these 
models, the RPM, and consequently the centrifugal force, was 
varied in the range (130 rad/s to 300 rad/s). 

The model geometric characteristics are the same as in 
Hodges[1987], a graphite epoxy composite rotor blade 


constructed out of Hercules IM6 fiber with Ciba-Geigy R6376 
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resin. The cured ply thickness is assumed to be 0.0055 in. The 
D=spar has a 35.23 in. radius, and@begins at station 532 in., 
which has all degrees of freedom suppressed to emulate a rigid 
rotor connection to the hub. The rotor has a constant cross- 
section, defined by 13 nodes and is divided into 21 span-wise 
(3.0 in apart) stations. Models with finer meshes were used 
in the early stages to check for convergence. The mesh 
selected seemed to be a reasonable compromise between the 
computational time and accuracy. 

The orthotropic material properties of the lamina are 


listed below: 


i; psi.....a... cet eee Doel sae? 
Bap) PSL ect ss sctseves 5 60 oh 6 aoe wed x mr0® 
Vigtvevvers es ; Py, NO iS 
ios DSi s+ +s sees Dear oece Oruececne 0.73 sae 0° 
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TABLE 1 = MODEL DERTNITION. 





Example GRIDS, POINTS AND LINES 
Vector 5 L45 4 L14 1 ou ie 
+20 
-70 
+20 iol 2 
-70 
-70 a2 Jaa eee 3) L8a 1 8 L128 2 
+20 (not to scale, e.g. line L512 is 30.0 in length, 
L45 is .286 in, and L14 is about 1.15 im) 
* = vector inward ; o = vector outward. 
POSSIBLE RIGHT HAND RULE NORMAL VECTOR TO SURFACE 
MODELS GRID DEFINITION 


GRID L512,L1112, | batt bse L18,L128,Lo129 
DEF. L411,L45// |L14// 1D de, 


GRID L512,b1112,) 047s eon Lo12,L128;, Digee 
DEF. ane, L4 5767) ba ody 


GRiv Lo12,L1112, | ble) boi, Laie L512,L128 , ee 
DEF. L411,L45// |L14// oe 


Rn ES TE EE ——— Le —————— —- ——— —— 


GRID L411,L1112,{/L18,L811,L411, L512,L128, biae 
DEF. L512,L45// |L14// Lode, 


Se er ee pn eee ee eee re — a eee - — 


a 


GRID L411,L1112, {b18, 081 ee L18,L128 Boies 
DEF. L512,L45// |L14// L51// 


ew mm ee ee 


GRID L411,L1112, | L411), 681 eee L18 ,L128, boas 
DETs. L512.) 7 


38 





ey STATIEC AND DYNAMIC 


ANALYSIS 


The six models described in chapter III are basically 
subjected to three main types of independent static loads. 
Further, a modal analvsis is also performed to obtain the 
first four fundamental modes of vibration frequencies and mode 
shapes. 

The first load applied is a torque at the tip, not 
combined with anv other load. 

The second load case consists of combined lift and blade 
weight at zero RPM, 1.e. with no centrifugal forces. 

The third load case is the centrifugal force due to 
rotation of the blade. The lift and blade weight are present 


and the RPM varies from 130 rad/s to 300 rad/s. 


A. TORS TONAL LOAD 


A moment of 26.4 in-lb is applied distributed among the 
nodes at the tip as concentrated moments. 

This load case reveals the influence of the laminate 
principal orientation angle (LPO), and the torsional stiffness 
for each model. The different response of the six different 


models also is presented for comparison purposes. 


Sus 


The results are presented in Figure 15 for the six models 
for LPO’s of +20° and -20° in each modem 
As much as 40% variation in the response may be observed 


in the comparisons. 
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Figure 15 Angular Displacement at the Tip Due to Pure 


Torque. 


| ee LIFT AND BLADE WETGHT LOADS 


This loading is designed to get the combined effect of the 


lift and blade weight given by Hodges[1987], expressed in 


terms of radial distance of the center as follows: 


(17) ) = (0.02222 lb/in) x r — OR 0R2 ae oan 
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This load is applied vertically in two spanwise rows at 
the leading and trailing edges, ina proportion that produces 
zero moments at the quarter chord. This load distribution is 
more homogeneous and induces less cross-section distortion, 
as opposed to the application along a line on upper surface 
as in beam type elements. This distributed line loading 
produces small torsional displacement (1074 degree). 

The results are presented in Figure 16 for the six models, 


with 2 values of LPO (+20° and -20°) for each model. 
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Figure 16 Displacement at the Tip Due to Vertical Loads 


Along the Blade. 
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ee. CENTRIFUGAL LOAD 


Next loading case investigated reveals the effect of 
varying centrifugal load, combined with changing the LPO angle 
within each model. The response vields twist angle induced by 
extension-twist coupling present in the structure due to 
asymmetric stacking sequence. 

The LPO, referred to the top layer angle, varies in the 
range of +/- 90°, the internal lay-up in the stacking sequence 
remaining constant, as previously defined. 

After some experiments to localize inflection points in 
the twist angle variation, the LPO to be investigated were 
chosen in the range of +/- 90°, giving a total of 11 laminate 


principal orientation angle for each of the six models. 
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Figure 17 Overall Range of the Twist Angle within the Six 


Models. 
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The loads due to lift and blade weight are present but 
their effects are observed to be negligible compared with 
those due to centrifugal loads. 

ine Seve salcCm ClmmesOommiacd/Ss to S00 rad/s in 
increments of 10 rad/s, in a total of i9 load cases. The 
combination of load cases and LPO give a total of 201 points 
of twiect angle for each of the s1x models. These points are 
presented in six 3-D graphs to better appreciate the choices 
offered in this desiaqn approach. Those graphs present the 
trends of the twist angle with RPM and LPO variations. The 
twist angle range for each model is in Fiqure 17. while the 
averall parametric information are presented in Figures 18 
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Figure 18 Model One, Twist Angle Induced by Rotational 


Inertia Loads Combined With Laminate Principal Orientation 


(LPO), Referred to Top Laver Angle. 


44 


MOvE 2 ee 


Sa ee 


wees. 
SRE © 


a a es SSR, 
ESOS 
i ee ea SO NS KF 


\ | 
en re Hh 
\\ 


tae 
> 
2% 
WS 
Ly 
® 
= 
Nt 
K 
eC 
Rs 
N 





Figure 19 Model Two, Twist Angle Induced by Rotational 


Inertia Loads Combined With Laminate Principal Orientation 


(LPO), Referred to Top Laver Angle. 
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Figure 20 Model Three, Twist Angle Induced by Rotational 


Inertia Loads Combined With Laminate Principal Orientation 


(LPO), Referred to Top Laver Angle. 
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Figure 21 Model Four, Twist Angle Induced by Rotational 
Inertia Loads Combined With Laminate Principal Orientation 


(LPO), Referred to Top Layer Angle. 
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Figure 22 Model Five, Twist Angle Induced by Rotational 


Inertia Loads Combined With Laminate Principal Orientation 


(LPO), Referred to Top Laver Angle. 
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Figure 23 Model Six, Twist Angle Induced by Rotational 


Inertia Loads Combined With Laminate Principal Orientation 


(LPO), Referred to Top Laver Angle. 
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JB es MODAL ANALYSIS 


To observe the effects of each configuration (Model) and 
laminate principal orientation on the first four fundamental 
frequencies, a free vibration analysis is performed. 

The results of the modal® analysis are depicted in Figures 
24 and 25 for the six models and for LPO’s of -20° and +20° 


respectively. 
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Figure 24 Frequency Variation among Different Models. 


> 1 FW means First Flapwise; 1 CW means First Chordwise 
and So on. 
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It mav be noted that, even for simple models, this 
solution procedure is a time consuming task, only recommended 
for final or quasi final configurations in the design process, 
or uSing a very crude mesh lay-out. 

It can be seen that, for each vibration mode, the 
frequency remains practically constant. This reveals a low 
sensitivity of the frequency to changes in model 


configurations, considering the first four fundamental modes 


and same LPO’s. 
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Figure 25 Frequency Variation among Different Models. 
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Vv. ANALYSIS OF RESUVUET 


In this chapter the results obtained for each model will 
be presented and several observations are inferred. 

The present research was conducted to get an insight of 
the trends implicated into the two questions posed in Chapter 
III, sections B and C. This study also provides some insight 
into answering the questions posed in sections B and C of 


Chapter memn-s 


What is the magnitude of twist available from a particular 


geometry and material system? 


What are the possibilities to model different structural 
configurations starting with the same number of layers, 
same lay-up design, same thickness and material 


properties? 


The first question is a design approach to structural 


analvsis. That approach will be elaborated in this chapter. 


Se 


The second question opened a vast field for research 
activity that will be addressed in next Chapter. 

It can be seen in Figure 14 that the models 1 through 6 
are obtained simply by changing one of the normal vectors 
associated with the laminate, (in one of the three faces of 
the spar). The effects of these changes in the structural 
behavior of the blade may be divided into two main categories, 


static and dynamic behavior. 


Ae DYNAMIC BEHAVIOR 


The free vibration characteristics are shown in Figures 
24 and 25. Figure 24 shows the frequencies for LPO of -20° and 
Figure 25 presents for LPO of 20°. The first four fundamental 
modes of vibration have very little sensitivity to the model 
type, and may be observed that the frequencies do not 
change, for the same or symmetric LPO (positive/negative). 

This behavior can be considered benign, or in other words 
auspicious, in a connotation that, once the designer arrived 
at the desired frequency values for the RPM range considered, 
the LPO or the model can be changed until a desired twist- 
extension coupling is achieved. 

These iterations can be done without looking very close 
to resonance problems, which must be checked anyway when the 


Pena! model is to be selected. 


Cy 
a) 


This assumption is made based on the validity of 
Yntema[1955] approach for these models. It is possible that 
the change in the stiffness matrix due to RPM affects the 
rotor built of laminated composite materials in a different 
way than isotropic rotors. Further research needs to be done 
along this direction to obtain comparable results for 
composite rotors. 

In Figures 24 or 25, each column corresponds to one mode 
of vibration at zero RPM and is a point in the 3=-D plotete. 
that specific mode (flatwise, chordwise or pitchwise’). The 
RPM and LPO may be varied as in the static case yielding 
carpet plots. The resulting plots would validate or not the 


assumption of Yntema’s values for composite blades. 


B. STATIC BEHAVIOR 


The static behavior is presented in Figures 18 through 23. 
The parameters considered are RPM and LPO as inputs and the 
twist-extension coupling manifested through the twist angle 
at the tip, as the output. This angle will be referred as TTA 
(Lip twist angle) in the subsequent discussion. 

The TTA can be used as a measure of the twist angle ratio 
along the blade span. The twist angle is not a linear function 


of the radius but may be represented by a second degree 


9 Pitchwise mode is a higher mode in this case and is not 
considered in this report. 





Permxiomial, such that the TTA is a good parametric factor to 
infer the twist angle magnitude along the blade. 

As mentioned before, the six models are obtained by 
changing the orientation of the normal vector associated to 
the laminate, in the three shells that constitute the spar, 
one at a time as shown in Figure 14. 

In that figure it can be verified that Model 1 and Model 
4 can be built with only one unbroken piece of laminate, 
wrapped around a mold. Further, one model is just the inverse 
of the other, in other words Model 1 is the inside out of 
Mode] 4. 

mre Moaels 2 7and 39 can also be considered as some sort of 
inverse, ina broader sense; Similar reasoning may be extended 
memmodels 3 and 6. 

With the above considerations, the six models can be 
Semen Zed in pairs, which 1S very convenient from the point of 
View of getting the information in a simple and concise form, 
so that the designer can gquichiy evaluate the model, keeping 
in view the design requirements. 

The second important consequence of this inverse concept 
Meethat, Model 1 and Mode] + can supply a spectrum of twist 
angle requirements, for a given loading condition. 

The importance of this reside in the manufacturing 
Meee ss. in that both models may be built by filament wound 


Meelmigve. This process if less expensive and nore reliable 


as it eliminates the filament discontinuity inthe @sheumee 
junction, and requires less accurate qual iGyecormen oe 

The convention used in the modeling process has as a 
aftermath that positive angles have active aerodynamic 


stabilizing effect, acting as a propeller governor device. 


1. Models 1 and 4 
Figures 18 and 24 exhibit the behavior of these two 
models. 
The models, which were referred as inverse to each other, 
also show the response characteristics and behave as expected. 
This mav be noted by observing the displacement for a 


given RPM and Laver Principal Orientation Angle in Figures 


2. Modelss 2 ands 
Those models have the same inverse relationship in 
the TTA value. The main difference is that these two models 
are stiffer than the other four, as shown in Figure 17. Figure 


21% shows these two models drawn together. 


3. ‘Modelss3aand G6 
The inverse relationship also holds here and the 
stiffmess values are higher than that found in models 1 and 
+ and lower than models 2 and 5 (Figure 17). See Figitpemeas 


[On tins) compara s6Gm,. 
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Figure 26 Models 1 and 4 static behavior comparison, inputs 


are RPM and Laminate Principal Orientation (LPO) angle. 


Bur put is Tip Twist Angle (TTA). 


Figure 27 Models 2 and 5 static behavior comparison 


are RPM and Laminate Principal Orientation 
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Figure 28 Models 3 and 6 static behavior comparison, inputs 
are RPM and Laminate Principal Orientation (LPO) angle, 
BeEput is Tip Twist Angle (TTA). 


a, 


Vy ae CONCLUSIONS AND SCOPE 


FOR FUTURE RESEARCH 


A. CONCLUS TONS 


This report presents a detailed analysis of a typical 
rotor blade, the focus being the static behavior in the 
presence of extension-twist coupling due to asymmetric 
stacking sequence of laminas. 

Different extension-twist coupling effects are obtained, 
by changing the ply orientation with respect to the body axes, 
while keeping the stacking sequence, thickness and internal 
lay-up sequence as invariants. 

These effects are presented in terms of carpet plots. The 
extension-twist coupling is measured through the angular 
displacement at the tip of the blade as function of RPM and 
Laminate Principal Orientation ( LPO ) angle variations. 

Free vibration analysis is performed for two LPO’s of the 
six different models. The first four fundamental frequencies 
show small sensitivity to these specific structural changes. 

The static analvsis revealed a sort of "“antisymmetric™ 
behavior within the models, allowing them to be grouped in 
pairs. Such behavior can be of use in the design process, to 


trim structural response to prescribed loads. 
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A structural analyst should be very careful in dealing 
With composite laminated structures, especially if the analyst 
is experienced only in structural analysis of isotropic 
materials. It was noted that it is very easy to be mistaken, 
mainly if there are asymmetries present within the laminate. 
It is very important to know how the local coordinates systems 
are established in the finite element program being used. 
These axes will determine the relative position of the 
laminate with respect to the structure. If asymmetries are 
present, it is possible to get different structural response 
and have different structural configurations, all starting 


with the same basic laminate. 


B. SCOPE FOR FUTURE RESEARCH 


Future research may be beneficial into two main areas. 

The first one being to validate these models with 
experimental data, providing the necessary confidence in 
further utilization of the method envisaged. 

Once the validation is done, the method can be corrected 
and trimmed for operational use. 

Then a field that calls for attention 1s a dynamic 
characteristics mapping, that is, with the RPM and LPO as 
parameters obtain the fundamental frequencies as output. 

This would be a 3-D composite material extension of 


matcema s [1955] work for beams of isotropic material. 
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$ (COMPOSITE AREA =.086525 IN*2) 

$ (COMPOSITE VOLUME -=2.5957 IN‘*3) 

$ (TUNGSTEN AREA -.01358 IN‘*2 => DIA=.1315 IN) 

$ (MODEL FOCUSED LENGHT=30 IN; TOTAL LENGHT=35.23 IN) 

ORMAT ,53 


1/5.07E5,23.166,.338,.73E4, .004418 
pm], 1.4646,1.4164// 


¢ 

$ Poaeeeeewee TNITION [+207-707+20/- 70/7 -70/ +203 
$ fee Neos = .00OS5 IN/LAYER * 6 = .O33 IN 
COMETH 


ieee i,i,i,i,i 
0055, .01375,20 
,0055, .00825,-70 
0055, .00275, 20 
0055, -.00275,-70 
.0055,-.00825,-70 
Nem@ao, -.01375,20// 


¢ 

$ ELEMENTS 

= GallS MUST BE GENEXATED IN SAME ORIENTATION 

+ Peeetbloe ee VIEW PCINT NORMAL TO SURFACE OUTSIDE (OR INSIDE) 
+ Te secwmestaNolNG TH1S REFERENGE GENERATE THE GRIDS FOR 

$ ORTHOTROPIC MATERIALS (FILAMENT WOUNDED THINKING) 

$ 

$ oops eclPE OF GRID DETERMINES THE REFERENGE FOR MATANG 
5 

Gee bo/1,1// 

GRIDS4 


(eee) 6 LG11,L18,L14 
Greemeee, Li 78,L5i2,L51 
eer, lL iiiz,t411,Ll45// 
END 


ANNEX 2 


FILE TO BE USED IN BUILDING THE SEVERAL 
EACH M$**.SRC FILES 


$ = LAYUP NUMBER 

** = ANGLE OF TOP LAYERSWI SESS seas 
1- LAYUP GRID DEFINI Pan —-- 6. ete eS 
ewes ole 

GRID4 


G1/L411,L811,018,L14 
G2/L1e.lize.tsi2-ts 
GS/LS12.U1112. eae eee 


M2**.SRC 

GRID4 

Gi/Lai stare se ed 
G2/L512,.128,L!8,L51 
ea/lsiei11i2 oe L45// 
M3**,SRC 

GI 1 [ey 

Gi/l18 vis es 


G2/ULSl 202.1267 Pane 
G3/L512.L1112,L411,L45// 


MO? SRO 

GRID4 

(57 Pog ee eels. ies 
G2Z/Loeeyel ze. Liew! 
ee) g lealel see eer lene eee 
mate SAC 

GRIDG 

Gi/lLis ke te l4it eis 


G2Z/LIB, 2226) Bae od 
62/4 11 oie a ea 


Wie we 5 Sisi& 

GRiDG 

gay 2 We Ae eae eer co gh des ee cee a ear 

B27 tee eee See 
63/41 lois ee ee Ie 


LA'Y,U ©’ elke 


LAYUPS IN 


NORMAL VETOR DIRECTION 





fem LAYER ANGLE -- em | EMSS WS 

HEHEHE HHH HER HEHEHE HEH HHH HHH HEHEHE HEH HH HHH HHH 
meg. SRC 

meme, .013/79,0 

.0055, .00825,-90 

.0055, .00275,0 

Gees — .002/5,-90 

meres, - .00825,-90 

Meee, -.01375,0// 


KEK HHH HH EHH HH HEHEHE HHH HEHEHE EH HHH EHH HH HHH 


Peo. OKC 

meee, -01375,10 
meee, .00U825,-80 
feoss,.c 00275,10 
Meee .O00275,-80 
Bees, .00825,-80 
meee .013725,10// 


5 ee ee ee eee eee ee 
Mae. or 

fee, .01 375, 20 

meee, .O0875,-70 

fee UO. 75,20 

a BOO 27s. —- 70) 

eee ,O0025,- 70 

eG. ., — mola /S ee 


ee ee ee ee eee ee ee ee ee eee 
| 

fee, 0 3/5, 20 
Gee , VOL 25. -60 
meee, . 00775 ,30 
Sete , 0G 275, -60 
Reeeeo - .008725,-60 
eee, - .01375,30// 


HHH HHH HHH RHR RHEE HHH HHH HHH HH HHH EHR HHH HH HH 
MS4O.SRC 

ggeee, .01375,40 

Geeear, .00825,-50 

Gee. O02 75,40 

Geer 00275, -Su 

fee .O0825,-50 

eee ,01375,40// 


=. ] 


Crt 


fs eO. Ske 
BUOSD = a. 7a ane 
OOS cI 2 OCSe a 6440) 
OOZ Sas 0) 
S0G2 7a 5-46 


BOOS or, 
2OCSS 


Baa 


TOSS a= 


MOC sai. 


OCS Zone 40 
~OLS 75,9077 


MMM BREE HH HHH HH HHH HH HEHEHE EHEHEHHHAHRHH HH 


M$60.SRC 
.0055, .01375,60 
0055, .00825,-30 
00275, 60 
,00275,-30 


mOCSs 
KOS Si 


TOOS oa = 
Ee Clare Wigs 


OO8 2 2 ae 
POS aye Oly 


5 ie ed ede ee ee ee ee ee ee ee ee ee ee i ee ee a ee ee ee ee ae eee ee eee 


M70. See 
,0055,.01375, 
00825, 
.00275,70 


~OUS ae 
,OG Sai, 
~OODen, 
2 Cheha a 
Oe Ser. 


KRAMER HH EHH 


—" “ s 


cee OL ts 
OS a ae 


—_ 


See 7a 


MEeC. SRC 


70 
=2 0 
an ay 
a a - 
a Os 
dem ‘o 


PG 


> ee de ee eee ee ee eee eee ee ie eee ae ee ee eee ee ee 


OOS Seo 2 Sy aoe 
SOO Gon = VOR 2 on see 

~ VOSS. <OOZ7 >a 7a8 
~OOSar 4 OO? Warne? 
OOSS > , OOS Zara 
eOODD, =. O1TS7 a see 
5 ie ee ee ee i a ee ee oe ee ee ee ee ee ee ee eee oe eo 
Mie 2 Son 
OOOO 6 . 0 Same) 
SOOaDo.. OOS a.e 

OO Soe. OO 27 ore 
BOOS aie =O Gai ae 
Cs err e Es Sia tas. 

5 ONG tel en Oe amen Ty, 


Me-S0. 
a) Che bx 

CA on fees (ee 
s wes? we 4 


te an q + ra 


’ 
i §] CN cS 


a i 
‘eh 6 ah 


RMR RH HHH RH HRM RH RHE RH RH HH 





me-B0.SRC 
0055, .01375,-80 
0055, .00825.106 
fas), .00275, -80 
Mee5>, -.00275,10 
,0055,-.00825,10 
Wess, -.01375,-80// 


ee Te Te EE TE EE SEE HE IE EE EE EE EE EEE EE EE EE 


fee7o.SRC 
Ne@e>, .01375,-70 
0055, .00825,20 
0055, .00275,-70 
.0055,-.00275,20 
.0055,-.00825,20 
Meese. -.01375,-70// 


Ht Fe Te HE EH Ee FEE HE EE EEE HE EE EE EE KE ER EE 


M#-60.SRC 
Meos>, .01375,-60 
.0055,.00825,30 

Nemes, .00275,-460 
0055, -.00275, 30 
Sess. -. 00825, 30 


“oo & a ~~ f=. Suimle saa: 
owt) fe OL 2/5, HOS / 


HERE ERE REE EERE EE DREHER ERE REPRE HEHEHE ED HHS 
Mf-S0.SRC 

Meee, .01375,-S0 

Geeee, . 00822, 40 

0055, .00275,-50 

Gee, - .900275,40 

Bees, .00625,4C 

Mees) - .01375,-50// 


2d ee ie ie i ie de ee ee ee ee ee an a oh cn ap a op aia o 
eee 9. SKC 

gem, .01375,-40 

fees, .00825,50 

meee, . O07 75,-40 

eee ~ .00275,50 

Gees - .00875,50 

wee, -.01372,-40// 


MT OE Oe te te te Oe ee eH ee 
Ma 50.SRC 

See, O1575,-30 

eee DOB LS, 60 


ie ani eee aa l= SNe th 
° US w~—- 5 8 woo 7 =. 4 a ae! 
in Ao es Lame ae Ss = 
eC ete. O Tae ie 7S, AG ape) 


—_ ~ 
— 


Oo S } = WORE Ay 


= — @ e ‘= er, 


fi) f 
ee, 


ee 6° SIO 


a ese ihe Sma cies mers § 


~~] 


KREME HERE RHEEREEHHHEHE HEHEHE 


M$-20.SRC 

0055, .01375,-20 
0055, .00825,70 
0055, .00275,-20 
.0055,-.00275,70 
.0055,-.00825,70 
10055, 861375, -26// 


sie ie ie ie ie ee ee ee Ee ee i ee ee ee eee Ee ee ee ee 


MSO. Sve 


.0055,.01375,-10 
0055, .00825,80 
,0055,.00275, - 10 
.0055,-.00275,80 
.0055,-.00825,80 
AOO@sa = O1s 754-1 oy, 


“wer eAb Pisa t Ra Burton 


Defense Technical Information Center 
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Library, Code 0142 
Naval Postgraduate School 
Monterey, CA 93943-5002 
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IPD - Divisao de Aeronautica 
Attn: Chefe da PAR 

Sao José dos Campos, SP BRAZIL 
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Department of Aeronautics 
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ieee 1 


N° Copies 
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ii 
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